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Abstract 

We explore the physics potential of multi-megaton scale ice or water Cherenkov detectors with 
low (~ 1 GeV) threshold. Using some proposed characteristics of the PINGU detector setup we 
compute the distributions of events versus neutrino energy Ei, and zenith angle 6^ , and study their 
dependence on yet unknown neutrino parameters. The {E^, — 6^) regions are identified where the 
distributions have the highest sensitivity to the neutrino mass hierarchy, to the deviation of the 
2-3 mixing from the maximal one and to the CP-phase. We evaluate significance of the measure- 
ments of the neutrino parameters and explore dependence of this significance on the accuracy of 
reconstruction of the neutrino energy and direction. The effect of degeneracy of the parameters 
on the sensitivities is also discussed. We estimate the characteristics of future detectors (energy 
and angle resolution, volume, etc.) required for establishing the neutrino mass hierarchy with high 
confidence level. We find that the hierarchy can be identified at 3a - 10a level (depending on the 
reconstruction accuracies) after 5 years of PINGU operation. 
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I. INTRODUCTION 



Atmospheric neutrino studies have enormous physics potential which has not been fully 
explored yet. Atmospheric neutrinos can be used (i) to explore different effects of neutrino 
propagation: oscillations in vacuum and in matter, MSW resonance of neutrino oscillations 
[HE], as well as parametric enhancement effects [3], etc.; (ii) to determine neutrino oscillation 
parameters: mass squared differences, mixing angles, CP violation effects, mass hierarchy; 
(iii) to search for new physics beyond the standard framework with 3 light neutrinos: new 
neutrino states, non-standards interactions, violation of fundamental symmetries, etc.; (iv) 
to perform, in principle, the oscillation tomography of the Earth. 

Among major results obtained with the atmospheric neutrinos are the discovery of neu- 
trino oscillations [4j, first measurements of 623 and Am|2 [S] and various bounds on new 
physics [6]. 

The physics potential of atmospheric neutrino studies with existing and forthcoming 
experiments has been widely discussed before (see, e.g., [7j and references therein). In 
particular, the potential of the IceCube's DeepCore (the currently existing detector within 
the inner core of IceCube with an energy threshold Ei, ~ 10 GeV [8J) for studying neutrino 
oscillations has been explored in Refs. P4TT]. The first experimental results on atmospheric 
neutrino oscillations in DeepCore have been reported in Ref. [T2] . 

Although at the probability level the effects of the neutrino mass hierarchy, deviation 
of the 2-3 mixing from the maximal one, and CP-violation in specific oscillation channels 
can be of order 1, there are a number of factors which substantially reduce the effects 
at the level of observable events. As a result, determination of the neutrino parameters 
becomes difficult if possible at all. The main challenges include (1) relatively low statistics, 
especially at energies above a few GeV; (2) the presence of both neutrinos and antineutrinos 
in the original neutrino flux and the difficulty of experimental separation of the signals from 
neutrinos and antineutrinos, especially in large (megaton scale) detectors; (3) the existence 
of both z/e and in the original fluxes; (4) a significant smearing of the signal over the 
energies and zenith angles, related to large uncertainties in reconstruction of the neutrino 
energy and direction. 

Recently, the idea has been entertained that large statistics of events that can be col- 
lected over a wide energy range in multi-megaton detectors with low (a few GeV) thresh- 
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olds, supplemented by relatively mild technological improvements, can alleviate or remove 
the above-mentioned shortcomings [T3HI3]. High statistics will allow one to make specific 
selections of events with different geometries from certain ranges of energies and zenith an- 
gles in order to enhance the sensitivity of the experiment to various neutrino parameters, 
to resolve degeneracies between those parameters, and to reduce the effects of uncertainties 
of the input parameters (such as the neutrino fluxes and experimental resolutions). Re- 
cent considerations and discussions of the multi-megaton ice detector PINGU (the proposed 
upgrade of the IceCube detector) (TU] show that this hope may actually be reahzed. 

In this paper we explore the possibilities of studying the neutrino parameters in multi- 
megaton scale ice or water Cherenkov detectors with energy thresholds as low as a few GeV. 
We study the energy and zenith angle distributions of events and their dependence on the 
neutrino mass hierarchy (ordering), the deviation of the 2-3 mixing from the maximal one 
and on the CP-phase. We identify the geometry of the events and the kinematic regions in 
the {E,^ — cos 9z) plane where the dependence on a speciflc neutrino parameter is maximal. 
We compute significance of measurements of these parameters and explore dependence of the 
significance on the energy and zenith angle resolutions of the detector, i.e. on the accuracy 
of reconstruction of the neutrino energy and direction. 

In our calculations we use some provisional characteristics of the PINGU detector, in 
particular the effective volume and its energy dependence. At the same time, we perform 
our analysis with an effort to make it as much as possible independent of specific experimental 
features, which are yet to be determined. 

The paper is organized as follows. In sec. II we summarize relevant information about 
the oscillation probabilities. In Sec. Ill we present distributions of events in the E,y — cos 6z 
plane as well as significance plots for determination of various neutrino parameters. In 
sec. IV we perform smearing of the event densities using neutrino energy and direction 
reconstruction functions. In this section we also study dependence of the significance of the 
determination of the neutrino mass hierarchy on the energy and angular resolutions (the 
widths of the reconstruction functions) of the detector. In addition, we discuss here the 
issue of degeneracy of neutrino parameters. Sec. V contains discussion of our results and 
conclusions. 
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II. OSCILLATION PROBABILITIES 



In this section we discuss dependence of the oscillation probabilities on the neutrino 
mass hierarchy and on the deviation of the 2-3 mixing from the maximal one. We consider 
the evolution of three neutrino flavors Uf = {ue, v^., VtY^ in the propagation basis, Vp^op = 
{ue, z/2, i^s)^, defined as Uf = U2zhJ-'prop- We use the standard parameterization of the PMNS 
mixing matrix, Upmns = U23ISU13U12, where Uij is the matrix of rotations in the ij— plane 
and Is = diag{l, 1, e*'') is the matrix with CP-violating phase. 

We consider the neutrino energy range E^, > (2—3) GeV which includes the 1-3 resonances 
and parametric enhancement ridges, and where the sensitivity to the neutrino mass hierarchy 
is expected to be maximal. Also in this range one expects small CP-violation effects, so that 
their degeneracy with the effects of the mass hierarchy is small. In this energy range we 
can neglect the effects of 1-2 mixing and mass splitting in the first approximation. Then 
the oscillation probabilities can be written in the following forms, which are convenient for 
discussion of the neutrino mass hierarchy [T7]: 

Pee = I- Pa, (1) 
-f/tte — Pe^ — ^23-Pa ; (2) 
Per = C23P4 ; (3) 

Pmm = 1 - ^ sin" 2^23 - sIsPa + ^ sin^ 2^23 v^l - Pa cos 0x , (4) 
1 1 



Pf,r = 2 sin 26^23 - S23C23PA - 2 sin 26^23 V 1 - Pa cos (px, (5) 
where Pa = l^esP ^nd 

(px = arg[A22Al-^]. 

Here Aij are the amplitudes of z/j — )■ uj transitions between the states of the propagation 
basis. From these formulas one immediately sees correlations between different probabilities. 
The amplitudes can be represented as 



^33 ~ l^ssk**^^^ ~ \/^ ~ Pac^'^^^ , ^22 ~ ^ 



«<P22 



SO that (px = 022 - 033- 

In the case of neutrino propagation in matter of constant density we obtain explicitly 

km. V + A 



'■>x 



( <P31 \ V + A 

arctan cos2^?J tan — H x, 

V 2 / 2 



where 

cos 2^13 A - V 



v/(cos 2^i3A - Vy + A2 sin^ 2^13 ' 



cos2e™ 



= a;-^(cos2ei3A - y)2 + A2 sin' 2^13 

and 



2E. 



For matter of constant density 



In the hmit ^^13 — t- 0, we obtain 



sin2 2^r3sm'^. 



^32 



At the MSW resonance, tan2 6'5|^ = 1, and the phase 0x is again given by ^ (up to a 
factor cos^ ^13 ~ 1). The phase 0x approximately equals 032 also at energies far above the 
resonance energy and coincides with ^32 in vacuum. Eq. ([6]) is therefore violated only in the 
resonance region, except at the resonance point itself. This is a consequence of the smallness 
of the mixing angle 9i^. 

The probabilities for antineutrino channels can be obtained from the expressions in eqs. 
1|-[5|) by substituting P4 — )■ P4 and 4>x 4>x, where 

Pa = Pa{V ^ -V), = MV ^ -V). 

In the case of normal mass hierarchy (NH), the 1-3 antineutrino mixing in matter is sup- 
pressed, consequently. Pa ~ and (px nearly equals the vacuum phase difference ([6|. 

In the approximation of zero 1-2 splitting, for the inverted neutrino mass hierarchy (IH) 
we obtain 

pIH _ pNH ,IH _ _7NH 

— 1 — 'rx 1 

where ~ 'P^^iY ~^ (see the Appendix). Therefore 

tdIH _ pNH pIH _ tdNH 

The equalities in eq. ([T]) receive corrections from non-zero 1-2 mass splitting and mixing. 

The information about the neutrino mass hierarchy is encoded mainly in Pa and also in 
(px- If -Pa = and (px = 032, the oscillation probabilities for the IH and NH coincide. 
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FIG. 1: Dependence of the oscillation probabilities in various channels on neutrino energy for a 
number of values of the zenith angle. Normal mass hierarchy is assumed. Left panel is for the 
Vx channels, and right panel is for the Ve — >■ Vx channels, where Vx = Ve (blue lines), (red 
lines) and (black lines) . We use the values of the neutrino parameters defined in the text and 
(5 = 0. 

Let us consider in more detail eq. (|4]) for the probability P^^, which plays a crucial role 
in our analysis (a similar analysis can be performed for P^,- as well.) The first two terms 
in (|4| correspond to P, the average 2z/ probability with vacuum oscillation depth, which 
is due to 2-3 mixing. The probability P^ is an oscillatory function of the neutrino energy 
(and the zenith angle), but in the resonance region the period of oscillations is determined 
by the oscillation length in matter ZJ^ l^zj sin 2^13 ^ (3 — 4)/23, which is much larger than 
the oscillation length Z23 responsible for 0x- The effects of the 1-3 mixing (and therefore 
the matter effects) (i) reduce the average probability, the third term in Q; (ii) reduce the 



depths of oscillations by the factor \/\ — Pa, the fourth term in (|4|); and (iii) change the 
oscillation phase 032 — )■ ipx- The third and fourth terms lead to a change of the depths of 
the oscillations. These features are well seen in Fig. [T| where we show the dependence of 
the oscillation probabilities in various neutrino channels on the neutrino energy for different 
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FIG. 2: Probabilities of oscillations of neutrinos of various flavors (blue lines for v^, red lines for 
z^e) to vs. neutrino energy for a number of values of the zenith angle. The solid and dashed lines 
are for NH and IH, respectively. Left panel is for neutrinos, and right panel is for antineutrinos. 
All the parameters are the same as in Fig. [TJ 

values of the zenith angle. Figure |2] shows the oscillation probabilities for both the neutrino 
and antineutrino channels and for the two mass hierarchies. 

We use the values of the neutrino parameters = 2.35 ■ 10^^ eV^, Am'^i = 7.6 ■ 10~^ 

eV , sin ^23 = 0.42, sin^ = 0.312 and sin^ ^13 = 0.025, which are close to the current best 
fit values [IH], and the PREM density profile of the Earth [12] for numerical computations. 

The strongest modifications of the oscillation probabilities due to matter effects are in 
the resonance region E,, {A - 8) GeV, where Pa has a peak due to the MSW resonance in 
the mantle of the Earth. For the selected values of the oscillation parameters the maximum 
of the peak is at E^, = 6.2 GeV and cos^^^ = —0.68. In addition, at this point the oscillation 
phase 0^ = 71. In the region | cos^^l > 0.83 and E^ < 1 GeV strong modifications of the 
oscillation probabilities are due to the parametric enhancement of the oscillations for the 
neutrino trajectories crossing the Earth's core, and also due to the MSW resonance in the 
core [131 [H]. Notice that P^^, and consequently the mass hierarchy effects, strongly depend 
on the 2-3 mixing: the probability Pa enters P^^ with the factor S23. 
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FIG. 3: Neutrino oscillograms of the Earth (lines of equal probabilities in the Ey — cosOz plane) 
for different oscillation channels and for the values of the oscillation parameters indicated in the 
text. Shown are the oscillation probabilities normalized by their maximal values in the parameter 
space of the panels: Pa^ / P^ff' ■, with P^""^ = PJ^^^ = 1. is in GeV. Normal mass hierarchy is 
assumed. 

In Fig. [3] we show the neutrino oscillograms - curves of equal oscillation probability in 
the {El, - cos 6z) plane - for the — )■ and z/g — ?■ Ux channels. The probabilities increase 
monotonically from darker-shaded to lighter-shaded regions. 

III. NUMBERS OF EVENTS IN THE NEUTRINO ENERGY AND ZENITH AN- 
GLE PLANE 

We consider first the numbers of events of different types produced by neutrinos with 
energies and zenith angles in small bins ^{E^) and ^{cosOz). These would correspond 
to real observables if the neutrino energy and zenith angle could be reconstructed from 
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the measurements with neghgible errors. In the next section we will study the effects of 
smearing of these distributions due to uncertainties of the neutrino energy and zenith angle 
reconstructions. 

A. Distributions of the f^-like events 

The z/^-like events (tracks) correspond to interactions v^ + N ^ ^^ + X, v^ + N ^ /i^+X. 
There are also some contributions from v-r which produce r with subsequent decay into /i. 
The number of z/^-like events in the ij-bin is 

Nf^^ = 2'kNapT I dcosO, I dE, V,^{E,)D^{E,,e,), (8) 

J Ai cos 6*2 JAjEi, 

where T is the exposure time, A^^ is the Avogadro's number, p is the density of ice, V^s 
is the effective volume of the detector, and the number density of events per unit time per 
target nucleon is given by 

D,{E^, e^) = [a^^ (<|.JP,, + $°Pe,) + a^^ ($JP,, + ^^P^,)] . (9) 

Here = ^'^{E^,9z), are the original fluxes of neutrinos z/^. We use the effective volume 
of PINGU with 20 strings [20] which can be parameterized as 

pVes{E,) = 14.6 X [log(E,/GeV)]^-^ Mt. (10) 

The volume increases from 2 Mt at -Ei, = 2 GeV to 20 Mt at E,y = 20 GeV. (In general 
Ves depends also on 9z.) We have found Pa and the probabilities Pa/s = Pai3{E^,9z) by 
performing numerical integration of the evolution equation for the complete 3z/— system. 
We use the deep inelastic cross-sections 

a^^{E^) = 7.30 ■ 10-^\E^/GeY) cm\ 
a^^{E^) = 3.77- 10"^^(E^/GeV) cm^ 

and we take the cross sections for electron and muon neutrinos of the same energy to be the 
same. In our calculations we use the Honda et al. atmospheric neutrino fluxes [21], which 
were calculated for the Kamioka site; however, for neutrino energies above a few GeV these 
should also give good approximations for the fluxes at the South Pole. 
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The fine-binned distribution of events ([8j) witli A(cos^2) = 0.025 and AE^, = 0.5 GeV is 
sfiown in Fig. |4j Tlie number of events decreases with Ei,. The pattern of the event number 
distribution follows the oscillatory picture due to the main — mode of the oscillations 
with a clear distortion in the resonance region. The maxima and minima are approximately 
along the lines of equal oscillation phases E^ ~ </)32Am32| cos9z\Re (where Re is the Earth's 
radius), again with a distortion in the resonance region E^ = A — GeV. In the high event 
density bins the numbers of event reach ~200, and the total number of events is about 7-10^. 

Introducing the ratios of the fluxes, 



$0 
$0' 



$0 

$0- 



we can rewrite the expression for the density of events (|9]) as 



where 



P + -P 



P + -P 



(11) 



Recall that the ratio r = depends both on the neutrino energy and zenith angle, e.g., 

in the range E^, = (2 — 25) GeV and for cos 6^ = —0.8 the ratio can be roughly parameterized 
as r = 1.2 ■ {E^/1 GeVf -^^. 



B. Hierarchy asymmetry 

Let us consider the effects of neutrino mass hierarchy on the distribution of the events. 
Using relations ([T]) we can write the density of the events for the inverted mass hierarchy in 
terms of the oscillation probabilities for the normal mass hierarchy as 

1 



D 



IH 



P 



-Pr 



(12) 



Then the difference of the numbers of events for the inverted and normal mass hierarchies 
equals 



rNH 



«j,M u.M 



27tNapT I dcose, dE,V,s{D'^-D 

'Ai cos 6*2 JA-jEi, 



where 



D 



IH 



D 



NH 



^1 l^fj) (-PjLt^t -^M/i) "I" ^ (1 ^e) {Pefj, Pefi) 



(13) 
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FIG. 4: The fine-binned distribution of the number of i/^— like events in the {Ei, — cosOz) plane 
that can be collected by the PINGU detector during 1 year; NH is assumed. 



and 



K 



e 

In the approximation Pa ~ 0, which is justified if the true hierarchy is the normal one, we 



obtain P^f^ ~ and 



P^ipi ^^~\ Si'^^ 26'23 [1 - cos 032] , 



where we have taken (px ~ 032 



. Consequently, 



P — P 



P — P 



- sin^ 2^23 
—sI^Pa- 



cos 



032 - a/1 - Pa cos 0x 



+ st.,P. 



23-^ A , 



(14) 



In this approximation 



1 - ^ sin^ 2^^23] (1 + - si, (si, -1]Pa+ 



sin^ 26*23 ( f^fi cos (^32 



1 — Pa cos ( 



>x 



The sensitivity of this quantity to the neutrino mass hierarchy is due to Pa 7^ 0, and so the 
highest sensitivity is expected in the kinematic region where Pa is relatively large. 
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For the difference of the numbers of events for the two hierarchies (13) we obtain 

^IH _ ^NH 



^^^"^-J I ^sin^ 2^^23(1 - ^^^) (cos(/)32 - ^I-Pacos, 



X + 



+ -^23 



,1 f^fi)^23 



f 



Pa 



Let us introduce the N-I hierarchy asymmetry for the ij-bin in the {E^, — cos plane as 



A 



N-I 



(15) 



The moduh of the asymmetries ( 15 ) are the measure of statistical significance of the differ- 
ence of the number of events for the normal and inverted mass hierarchies: Sij = \Aij\. 

Let us consider the condition Nj^^ = Nf^^ which gives the borders of the regions in the 
{El, — 9z) plane where the difference of the numbers of events has definite sign. It coincides 



approximately with the condition D 



IH 



D^^. The latter determines the lines of zero N-I 



hierarchy asymmetry. Using eq. (13) and approximate expressions in (14) we find from this 
condition 

cos (f)s2 - a/1 - Pa cos (f)x = - ■ s 



2^23 



r 1 - K,, 



23 



(16) 



The phases ^32 and (^x are functions of Ey and Q^- Since cos 032 varies with {cosOz/ Ey) much 
faster than r{Ey,6z), it is this periodic function that determines the lines of zero hierarchy 



asymmetry. Our calculations show that eq. ( 16 ) determines the zero asymmetry lines rather 
well. 

In Fig. |5] we show the values of the hierarchy asymmetry in the E^ — cos 6^ plane. The 
maximal asymmetry is achieved in the low energy parts of the domains bounded by the lines 



of zero asymmetry (16). The large asymmetry is in the strips along the constant phase lines 
in the energy interval E^ ^ {A — 12) GeV, where these lines are distorted by matter effects. 
The asymmetry changes the sign with changing zenith angle. For instance, in the energy 
range (7 - 11) GeV there are three distinct zenith angle intervals with the asymmetry sign 
being the same within each interval. The number of such intervals increases with decreasing 
energy. Therefore, due to the averaging, the region of high sensitivity to hierarchy will 
shift to higher energies if the zenith angle reconstruction becomes worse. Furthermore, the 
regions of high significance of the hierarchy determination overlap substantially with the 
regions of small numbers of events (see Fig. 111). This means that the significance is enhanced 
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FIG. 5: The hierarchy asymmetry of events in the E^, — cos 9z plane. The absolute value of the 
asymmetry in a given bin determines the statistical significance of the difference of the numbers of 
events for the inverted and normal mass hierarchies. 



due to the smallness of the denominator in eq. (15), and it would be diluted by combining 
a given bin with bins which have higher statistics but smaller significance. 

There is an important background to the i/^ events which comes from Ur interactions. 



C. Vr events 

The Ur flux appears at the detector due to z/^ — z/^ oscillations. In Fig. |6] we show the 
distribution of the Ur CC events + — )■ r + X) in the E^, — cosOz plane. The figure 
is a kind of inversion of Fig. |4], with maxima substituted by minima and vice versa. The 
number of events is, however, smaller than the number of z/^ events due to the smaller cross- 
section near the threshold. Notice that since u^- (as well as ui) from the sequential r decays 
(r -> z/^ + X, T — !■ z/,- + / + z//) is not detected, the energy of the original Ur cannot be 
reconstructed. 

The u-r interactions 

u^ + N-^T + h-^fi + u + u + h (17) 
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FIG. 6: Same as in Fig. |4j but for u-r CC interactions. 

will contribute to the main sample of z/^ events with a muon and a hadronic cascade in the 
final state. However, the number of these events is relatively small, and in addition these 
events have certain features which can be used to discriminate them from the true events. 



Indeed, on average the two neutrinos which appear in the process (17) will take about 



1/3 of the energy of the initial neutrino. Therefore, for a given observed total energy 



+ Eh, the energy -Eo('^r) of the original neutrino in the process (17) should be about 
1.5 times larger than the energy EqIu^) of the true z/^ event: Eq{i't-)/Eq{i'^) ^ 1.5. If 
only Ef^ is used to reconstruct the energy of the original neutrino the rescaling coefficient is 
Eo{vr) / Eq{v^) = 2.5. On average we can take a factor of 2 for our estimates. Furthermore, 
the branching ratio of tau decay into muon is 5^ = 0.17. Consequently, the number of z/^ 



events due to the reaction chain (17) is suppressed with respect to that of the true CC z/^ 



events of the same energy by a factor 

Hm) I '''' 

(provided that the initial z/^ and Vr fiuxes are equal). The power of the second factor follows 
from the energy dependences of the neutrino fiux (oc E~^), the cross-section (oc E) and the 
effective volume (which we take to be oc E^'"^ here; note that in the energy range E^ = 10 — 35 
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GeV such a simple power law approximates eq. (10) within 5% error). The last factor takes 
into account the threshold effect for the r production. For the energy rescaling factor of 1.5 
- 2 we obtain a suppression factor for the number of i/^-induced z/^ events to be 0.05 — 0.08. 
Due to the missing energy and momentum taken by the two neutrinos in the final state of 



reaction (17), the smearing effects in the energy and angle of the original neutrino around 
the average values will be stronger than for the true z/^ events. Therefore, the observed 
energy + will have a bigger spread for z/T--induced i/^ events. Furthermore, they will 
be characterized by a larger average angle between the momenta of the original neutrino 
and the muon. 



There are other properties of reaction (17) which can be utilized to disentangle it from 
the CC z/^ detection reaction. In particular, correlations between E^ and Eh are different for 
these two cases. Furthermore, one can select E^, — 6^ regions in which the z/^ — transition 
probability, and consequently the z/,- flux, are suppressed. Such regions can be readily found 
with the help of Fig. [6] The corresponding restriction of the E,^ — 9^ parameter space will, of 
course, result in a loss of the overall statistics, but would provide us with cleaner events. The 
resulting statistics loss should be affordable because of the extremely high overall statistics 
in multi-megaton detectors. 

In principle, one could also sum up the z/^ and Vr events and consider them in the 
6*^ — {Efj_ + Eh) plane. However, to determine whether or not a useful information can be 
extracted from such data would require an additional analysis which is outside the scope of 
the present paper. 

For the above reasons, in what follows we do not explicitly consider the contributions of 
the z/t- — )■ r — J- /i events and simply treat them as a 5% systematic error. 



D. Cascade events and the mass hierarchy 

Following the IceCube terminology, we will call the events in which the muon track is not 
identified as cascade events. There are several different contributions to the cascade events, 
including even events with faint muon tracks which can not be identified. All z/^-induced 
CC events, except those when the tau decays into a muon, contribute to the cascade events. 
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FIG. 7: Same as in Fig. |4]but for i/g CC interactions; NH is assumed. 

For the cascade events Uf. + N ^ e + X and Ue + N ^ e'^ + X we have 



2nNApT 



d cos 9z 



Aj cos I 



dE,Ves{E,)De{E,, COS 9,), 



where 



D,{E, COS^J = a^^<l>° [(Pee + rP^e) + /«e (^ee + rP^e)] • 

In terms of the probabihty Pa the number density of events can be written as 



23 



1)1}. 



The event distribution is shown in Fig. [7j Notice that here the distribution is weakly affected 
by the oscillations due to a substantial screening effect: the oscillatory part of the number 
of events contains terms proportional to rs^s — 1 and fs^s — 1 which are nearly zero at low 
energies. 

The difference of the number densities of events for the inverted and normal mass hier- 
archies is 



^IH _ ^NH ^ ^CC^O (p^ _ p^) [(,^2^ _ 1) _ ,^(,,2^ _ 

The expression in the square brackets here can be rewritten as 

(1 - Kp) = (1 - K 



^■523 



(19) 



(20) 
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Thus, there is a double suppression of the difference of numbers of events: (i) due to the 
neutrino-antineutrino factor (1 — k^) related to the presence of both the neutrino and an- 



tineutrino fluxes, and (ii) due to the flavor screening (the last term in eq. (20)) related 
to the presence of both Ue and z/^ in the the original atmospheric neutrino flux , more 
precisely, due to the ratio of these fluxes being close to 1/2 at low energies. The difference 
— can be further suppressed by the smallness of Pa (or of the difference of the 
neutrino and antineutrino probabilities Pa — Pa)- 

The numbers of events for IH and NH are equal, Nj^^ = N^^, in the bins for which both 



sides of eq. (20) vanish. From this we obtain 

sl^{r - fKe) = (1 - Ke). 



For r ^ r it gives 



■^23 



So, in this approximation we have only one line of zero asymmetry. 

The hierarchy asymmetry in the z/g CC events is shown in Fig. [8j Maximal asymmetry is 
in the resonance region E,^ = {A — 8) GeV where it has essentially the same sign, so that the 
suppression due to averaging is absent. Unfortunately, other contributions to the cascade 
events have different pattern in the E,^ — 6^ plane. The number of cascades induced by u^, 
quickly decreases with the increase of the neutrino energy. 

The dominant contribution to the cascade events comes from the Ur flux which appears 
at the detector due to the oscillations: 



^NH ^ ^C^C(^^)$0 I 1 gj^2 2^^^^^ _ ^l_p^cOS0x] - cl,PAis 



2 

23 



1/r) 



For the z/^ events the screening due to the flavor composition of the original neutrino flux is 
absent. 

The neutral current interactions of all neutrino species contribute to the total numbers 
of the cascade events but do not affect the NH-IH differences of events. This reduces the 
hierarchy asymmetry: 



cascades 



A ^ e+T e+r 

■ri, 



As we will see, analysing only events will be sufficient to establish the neutrino mass 
hierarchy. We therefore do not include cascade events in our discussion. Clearly, cascade 
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FIG. 8: Same as in Fig. [5] but for Ue induced events. 

events can give additional information, but taking them into account would require a more 
sophisticated analysis. Notice that possible identification of r events in a large liquid argon 
detector and its consequences have been discussed in [22]; should such an identification turn 
out to be possible also in PINGU, it would increase PINGU's sensitivity to the neutrino 
parameters. 



E. Effects of deviation of the 2-3 mixing from the maximal one 



We describe deviation of the 2-3 mixing from the maximal one by 



C?23 = J - S 



23- 



From eqs. (|lT|), g and Q we find 



D 



NH I 



'23) 



1 - Pa cos (j)xK.^{l - cos 032) 



+ d 



23 



^23 I 



(21) 



Notice that both terms in (21) are positive for 623 < vr/4 



In Fig. |9] we plot statistical significance of the determination of a deviation of the 2-3 
mixing from the maximal one. Here again high significance regions coincide with the regions 
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FIG. 9: Statistical significance of the determination of tlie deviation of the 2-3 mixing from the 
maximal one. The difference of the tracking events for two values of the 2-3 mixing: sin^ ^23 = 0.5 
{Nf^f) and 0.42 (iV^^^'')- 

of low density of events (compare Figs. [9]and|4]). This means that an integration over large 
— cos 6z domains would lead to a dilution of the significance. 

F. CP-violation effects 

We can write the oscillation probabilities in the presence of a Dirac-type CP-phase S as 

-Pq/3 = PaP + Pip ■ (22) 

Here and P^^ are the (5-independent and ^-dependent parts of the oscillation probability 
Pq^, respectively (note that 7^ P^^p^)- ^ matter with symmetric density profile one 
has Ppa = Pap{S -6). 

Now the 1-2 mass splitting and mixing should be taken into account. The state i>2 does 
not decouple and the oscillation probabilities depend on the matrix elements A^2j ^23 of the 
evolution matrix A in the propagation basis. Since the amplitude is doubly suppressed 
(by small Amgi/Amg^^ and S13), the terms that are quadratic in can be neglected. We 
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have then 1151: 



p6 



P. 



P. 



0, 

sin 26'23|^e2^e3l COS {(f) + 5) , 

— sin 26*23 cos 5 {174^2^4^21 coscp + cos26'23Re[74ig(y43g — ^422)]} 



(23) 
(24) 



where ch = aig(A*-A -). We will also use the notation \A -I 

^ "^^ ez e3' ' ez' 

now depends also on the parameters of the 1-2 sector). 



\A 



(where Pa 



Note that the last term in the curly brackets in eq. (24) is small if the 2-3 mixing is 



sufficiently close to the maximal one, and in addition the amplitude is small. We shall 
therefore use for our estimates the approximation 



cos2^23Re[Ai3(A.3 - A.,)] ^ . 



(25) 



The 5-dependent part of the expression for the number density of the /i-like events is 



p,.„ + -p: 



MM 



pS _pS 
MM ' ^ EM 



(26) 



Next, we notice that in the case of normal hierarchy one has P5 ~ 0, ~ 0, and in the 



approximation (25) we therefore have P^^ ~ 
find 



PL ^ 0, P^^ ^ 0. From eqs. (|23|) and (|24D we 



Di-Di='' =a^^$° sin 2^23 /fV^ 



r - 1 1 

cos 0(1 — COS 6) sin sin 6 



(27) 



This, in particular, means that the difference Nf^ — Nf^ ° should nearly vanish whenever 



Ps = or Pa = 0, i.e. along the so-called solar and atmospheric "magic" lines [151 123^ 
The vanishing of the difference A*"^ — A'^^^^ is, however, not exact, as it relies on the 



approximation (25). In Figs. 10 - 11 we show statistical significance of measurements of 
the CP-phase. The strongest effect is at low energies: ~ 3 — 5 GeV. Notice that with 
increasing 6 the size of the asymmetry increases, but the regions of different signs of the 



asymmetry do not change. This is in agreement with eq. (27). Indeed, in this equation the 



dependences of the right-hand side on and 6 effectively factorize, because in most of the 
parameter space either the first or the second term dominates. For the z/e-like events we 
obtain similarly 



Dl ^ [(Pi + rP;^,) + [Pi + fPl)] ^ a^^<l>° r cos(0 - 5) . 



(28) 
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COS 6^ 

FIG. 10: The difference of the numbers of events for 6 = it/2 and 6 = 0. Statistical significance 
of measurements of the CP-phase. 

Here we have taken into account that Pee and Pee are ^-independent and that P^e is strongly 
suppressed in matter. For the difference of the densities of z/g-hke events we then find 

_ Dl=° = a^^^^r sin 2^23 VPaPs [cos (/>(cos 5 - 1) + sin sin 5] . (29) 

The hues of zero Pa and Ps (the atmospheric and solar "magic" lines) form a grid, which 
leads to a domain structure in the {E^, — cosOz) plane. For maximal 2-3 mixing the domain 
structure of the distribution of the events becomes sharper (Fig. [12]). As discussed in ref. [T5] . 
in the full Su framework there is no crossing of the "magic" lines determined by the solar, 
atmospheric and an additional phase condition, and instead there is a smooth transition of 
lines of different types into each other. 

IV. DETERMINATION OF THE NEUTRINO MASS HIERARCHY 

The fine-binned event distributions computed in the previous section allow us to identify 
the regions of high sensitivity to the neutrino mass hierarchy as well as to other neutrino 
parameters. They show the lines of zero asymmetry which separate the regions of same 
sign asymmetry. These distributions also allow one to identify the regions of high and 
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FIG. 11: Same as in Fig. 10 but for the phases 6 = 7r/4 and 5 = 0. 



low degeneracy between different parameters. By using them one can select the regions 
of integration over the neutrino energy and zenith angle in such a way that the hierarchy 
asymmetry is enhanced while the effects of the uncertainties due to the incomplete knowledge 
of the other neutrino parameters are canceled or suppressed. Obviously, one should avoid 
integrating over regions with different signs of the hierarchy asymmetry. 

A. Summation of events over different bins 



Since one needs to integrate over certain regions of the neutrino energy and zenith angle, 
the significances of determination of the neutrino mass hierarchy and of other neutrino 
parameters will be modified in comparison to those found for small individual bins. The 
combined statistical significance resulting from a summation over n bins is given by 



If the individual N^^ do not differ significantly, we have approximately 



n 

In ^-^ 
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FIG. 12: Same as in Fig. 10 but for phases 6 = tt/2 and 6 = and sin 



723 



0.5. 



The asymmetry has opposite signs in domains separated by the hnes of zero asymmetry, 
and the ultimate sensitivity to the neutrino mass hierarchy can be estimated as 



NH 



One can first estimate independently the sensitivities of individual domains with the same 
sign of asymmetry and then sum over the domains. The real sensitivity will actually be 
lower because of (i) the integration (smearing) over bins with different numbers of events 
and significances (which will dilute the significance of the most significant bins) (ii) the 
integration (smearing) over parts of domains which have opposite sign of the asymmetry, 
(iii) uncertainties of the other oscillation parameters, (iv) degeneracy of parameters, (v) 
systematic errors, etc.. We address some of these issues below. 



B. i/^— like events and the mass hierarchy 

The z/^— events produced by the charged current z/^ interaction are observed as muon 
tracks accompanied by hadronic cascades. For these events the energy of the muon i?^ and 
the direction of its trajectory characterized by the angles 6^ and 0^ as well as the total energy 
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of the hadron cascade Eh can be measured. Using this information one can reconstruct the 
neutrino energy: 

E^Ki Ef, + Eh-niN , 

where mjv is the nucleon mass. In fact the calorimetric (time integrated) measurement in 
IceCube provides directly measurement of Ei,. 

The reconstruction of the neutrino direction is more comphcated. In the first approxi- 
mation (at sufficiently high energies) one can simply use 91 ~ 9^ with certain spread which 
depends on the neutrino energy. More precise determination is in principle possible if one 
uses also the information about the hadron cascades. Knowledge of the energy of the cascade 
narrows down the uncertainty in the neutrino direction. Further improvement would be pos- 
sible if one determines the plane in which the muon and the original quark were propagating. 
In this plane one can introduce the angle between the muon and neutrino trajectories, 9^, as 
well as the angle between the directions of the quark and neutrino momenta 9q. Then, using 
the energy and momentum conservation laws and excluding 9q, one obtains the expression 
for the reconstructed neutrino angle: 

El' + El-El 



cos 91, 



'^ElE, 



where we assume that the muons are ultra-relativistic. In turn, the knowledge of dv and 
of the muon angles and 0^ would allow one to reconstruct the neutrino zenith angle 
9l^9l{9^,9,,4^,). 

There is a number of uncertainties in this reconstruction procedure: (i) Errors in the 
measurements of E^, 9^ and 0^; (ii) the uncertainty in the point of the neutrino interaction 
(i.e. of the beginning of the muon trajectory); (iii) the uncertainty in the position of the 
center of the hadronic shower and (iv) the error in the determination of the energy of the 
hadron shower. 

We will describe the uncertainties of reconstruction of the neutrino parameters by distri- 
bution functions for the reconstructed neutrino energies and zenith angles: 

GE{El,Ey), Gg{9l,9z), 

where Ei, and 9z are the true energy and zenith angle of the neutrinos. The distributions 
are normalized in such a way that 



/ 



dy Gyiy^y) ^1, y^E^,9^, 
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where the integrations are performed within the appropriate ranges of the parameters. For 
Gy we adopt the Gaussian form 

N 

Gy{y.ay) = -=^e ^ 
y iTiay 

where A^j^ is the normahzation constant, and cr^ and ae are the widths of the energy and 
angular reconstruction functions, respectively. Both widths depend on the neutrino energy. 
In this way we obtain the unbinned distribution of events in the {E^ — cos 6*^) plane: 

D,(E^cos^O = j dcose^J dE, Ge{EI,E,) Ge{e:,e,) Kff N^{E,,cose,), (30) 

a = e, fi. Then the binned distributions of events are 

N^I = 2ttNaTp [ dcosei [ dElD^{E\cosei). (31) 

We will explore the dependence of our results on the widths aE and ae- For the ideal 



resolution, Gy{y'^ ,y) = 6{y^ — y), we get from (30) and (31) the same results as before. 

Consider the opposite limit of large widths, 2ay ^ Ay. In this case it is worthwhile to 
interchange the integrations dE^d cos 61 and dE^d cos 6 z- Then for the box-like distribution 



functions, Gy = l/2ay in the intervals y = y'^ ± ay, formulas (30) and (31) reproduce the 
results for large bins Ay ~ 2ay. 

For a contained event (both the vertex and fi track are contained within the detector) 
the error in the reconstructed neutrino energy scales linearly with energy, i.e. aE„ ~ xEj, 
below ~ 100 GeV p8j. The error in reconstructing the neutrino arrival direction at low 
energies is hmited from below by the root mean square value of the scattering angle, 9rms ~ 
y^mp/E^, which corresponds to 17.5° at 10 GeV. 



In Figs. 13 - 16 we show the hierarchy asymmetry in the distributions of the z/^ events 
smeared with energy-dependent Gaussian reconstruction functions characterized by different 
aE and ae- After smearing we integrated the event density over the reconstructed energy 
and zenith angle bins of the size A(£'^) = 1 GeV and A(cos^2) — 0.05. The smearing leads 
to a substantial decrease of the sensitivity to the neutrino mass hierarchy. This reduction 
is a consequence of the integration over regions with different significance and statistics as 
well as over the regions with different signs of the asymmetry. 

Considering the effect in each bin as an independent measurement, we can find the 
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combined significance as 



2 



\ 



4 



where the sum over bins can be substituted by the integraL For illustration we assume that 
the uncorrelated systematic errors are proportional to the number of events: acorr = fNjJ^, 
where / depends on the binning. We will use / = 5% and 10%. In general / is a function 
of neutrino energy and zenith angle. Therefore the total error in each bin is given by 



Notice that since here the contribution from the systematic error is proportional to A^^, for 
the same / the role of systematic error decreases with decreasing size of the bin. 

Correlated systematic errors, e.g., those of the overall flux normalization and of the tilt of 
the spectrum, apparently can not reproduce the profile of the distribution of events similar 
to the difference of distributions for the two hierarchies. Therefore, their effect to a large 
extent can be reduced to a reduction of the exposure time and statistics. Moreover, these 
correlated systematic errors can be parameterized and reduced with better measurements 
of the flux. 

Instead of 6,y, one could consider the angle which is measured directly. Smearing over 
the angle O,, with ag ~ ^^/n^JE^ essentially corresponds to the transition from 6,^ to 9^. 
Again, due to measurements of the cascade energy this uncertainty will be further reduced. 



In Fig. 13 we used aE = Q.2Ey and oq = ^Jnip/ E^. It can be seen from the figure that 
the region of the highest significance is between the lines E^/GeV ~ 20(| cos^^l — 0.2) and 
Ey/G€V ~ 15(1 cos^^l — 0.2), and between cosdz = —1 to —0.7. It is shifted towards higher 
energies compared to the un-smeared case in Fig. [5j 

Detection of hadronic cascades can in principle improve the determination of the neutrino 



angle. In Fig. 14 we therefore use the angular resolution ae = O-b^ynip/E,^. Notice that with 
better angular resolution the region of high sensitivity shifts to lower energy and shallower 
zenith angle bins, approximately along the lines of high significnace. 

To study the effect of the energy dependence of aE on the neutrino mass hierarchy 
determination, in Fig. [15] we use the fixed energy resolution aE = '2 GeV and the angular 



resolution ae = ^Jm^jEy. A comparison with Fig. 13 shows that fixed aE leads to a 
somewhat higher significance of the hierarchy determination. 
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FIG. 17: Significance of hierarcliy determination vs. exposure time for various smearing schemes 



illustrated in Figs. 13, 14, T5| and 16 Different curves correspond to different systematic uncer- 
tainties (/) assumed in addition to statistical uncertainties. 



In Fig. [16] we used the energy resolution cr^; = 2 GeV but the angular resolution ae 
^.h^JrUp/E^. Figs. 



16 



and 



14 



show that with improving angular resolution the significance 
of the neutrino mass hierarchy determination increases significantly. 

Figure 17 shows the significances of the hierarchy identification, S**"*, for different smear- 



ing schemes described above and for different uncorrelated systematics (/). The 5-year 
significances from the plots can be compared with the results found when no smearing 
is performed: 5*°* = 45.5(t (no systematics), 5*°* = 28.9(t (/ = 5%) and 5*°* = 18.8a 
(/ = 10%). Note that if IH is the true hierarchy, then in the first approximation the re- 



sults simply correspond to Figs. 13 - 16 with reversed signs of the asymmetry, although 



significances become somewhat lower in general. 



C. Effects of parameter degeneracy 

In our computations of the hierarchy asymmetry we used the fixed values of Amgg, ^13, 
^23 and 5. In certain kinematic regions, uncertainties in these parameters may lead to the 
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FIG. 18: Effect of uncertainty in Amg^^. Shown is the difference of numbers of events for Am^j^ 
and Am|i + la = 2.47 • 10"^ eV^. NH is assumed. 



same effects as the hierarchy change. That is, the difference of the event distributions for 
the true and assumed values of the parameters may have the same pattern as the difference 
of event distributions for the normal and inverted mass hierarchies. The simplest way to 
reduce the degeneracy is to select only parts of the E^, — cos 9^ space where the effect of 
hierarchy change dominates over other effects. 

1. As follows from the comparison of Figs. 10, H , 12 and Fig. [5} the effect of CP-phase 
6 at high energies is smaller than 10% of the hierarchy effect, and so it can be neglected 
in the first approximation. Also, it is characterized by a different pattern, and therefore 
one can select the regions in the — cos 9z plane in such a way as to further suppress the 
effect of the CP-phase (e.g. the selected bins should include domains with different signs of 
CP-asymmetry) . We computed the hierarchy asymmetry for 5 = 7r/2 and found that indeed 
the CP effect can be neglected in the first approximation. 

2. Effects of the uncertainty in 6*23 are strong in the region along the line Ey/GeV ~ 
23 1 cos ^2 1 (Fig. [9]). In the high energy part of this region, E' > 10 GeV, the hierarchy 
asymmetry is small. There is still an overlap of the regions of large effects of the hierarchy 
and ^23 uncertainty at low energies. The bins can be selected in such a way that the 
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effect of the uncertainty in 623 substantially cancels out. Furthermore, one expects some 
improvements in the determination of 623 from the accelerator and combination of reactor 
and accelerator experiments, so that the degeneracy will be further reduced. In particular, 
combining the reactor data with the forthcoming results from T2K and NOvA will allow us 
to determine sin^ 623 with an accuracy between about 0.04 (for maximal 2-3 mixing) and 
0.008 (for sin^ 623 = 0.4) This has to be compared with the current uncertainty of this 
parameter 5(sin^ ^^23) — 0.05 |27j (all numbers correspond to la). 



3. The effect of uncertainty in Amg^ is illustrated in Fig. 18 , where we show the difference 
of numbers of events for the best fit and shifted upwards by la values of Am'^^. (A downward 



shift by la switches signs in Fig. 18, although with somewhat smaller significances.) In the 
limit Am2i = variations of Am|]^ are equivalent to a corresponding shift of the oscillatory 
pattern in the energy scale (see Fig. [2]) at high energies {E > 8 GeV). Note that the shift 
is different for different zenith angles and for neutrinos and antineutrinos. The regions of 
the strongest effect of this shift on the hierarchy asymmetry have substantial overlap with 
the regions of strong hierarchy asymmetry, and therefore certain selection of the integration 
regions (binning) is required in order to reduce the degeneracy and disentangle the two 
effects. Clearly, further improvements of the accuracy of measurements of Aml^ by MINOS, 
T2K and NOvA will alleviate this problem. In particular, T2K and NOvA will measure 
Am'^i with an accuracy of 5 x 10~^ eV^ (la) [26]. This is about a factor of two better than 
the current uncertainty of this parameter 



To further explore if the uncertainty in Amg^ can mimic the "wrong" hierarchy, we 

to the best-fit 
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applied the same as = 0.2Eiy and ae = \/ nip/ By smearing as in Fig 
and +la deviation of the Am|^ values. The resulting significances are plotted in Fig. 19 
Note that the region of the highest sensitivity to Am|^ has shifted to higher energies and is 
concentrated in a narrower, cos 6'^ < —0.8, range as compared to the region of the highest 



significance of the hierarchy asymmetry (cf. Figs. 19 and 13). 



Figure |20] shows the difference between the significances plotted in Figs. [19]and[l3j Notice 
that the significances of the hierarchy determination in the region i?,^ ~ 6 — 13 GeV are 
reduced considerably. The effect of the degeneracy can be suppressed by making summation 
of significances only over certain domains in the Ey — cos 9^ plane. 

To estimate the effect of the uncertainty of the value of Am\^ on the significance of 
the hierarchy determination, we simulated the data for NH A^^''^ for a fixed "true" value 
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FIG. 19: Same as in Fig. 18 but after smearing the t'^ events in the {El-cos 9^) plane with 



aE = 0.2£'iy and ae = -^/ra^lEy. 



of Am|]^ j^^g and then fitted them with IH, treating Amgs as a free parameter. We then 
minimized 5**°* with respect to Am32 and found the corresponding values of Am23 and 
S'^*„. Next, we repeated the same procedure for different true values of Amg^^^^^g within its 
Icr allowed range. This procedure is illustrated in Fig. |2T| The left panel shows 5**°* versus 

2.35 ■ 10^^ eV^ (vertical line). In the right panel of Fig. 



Irue 

the values of S'^*„ obtained through this procedure as functions of Am 
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we show 



31, true 



(solid lines). 



For comparison we show also 5**°* for Am23 fu 



Am|]^ true' without variations of Am^ 



(dashed lines). As follows from the figure, variation of IS.m^^ reduces the significance of 



the hierarchy identification 5**"* by ~ 50%, and this reduction weakly depends on Amg^ 



true' 



The values ^fn^-i^jit^min within the present la uncertainties of determination of this mass 



difference (2.17 ■ 10 ^ — 2.59 ■ 10 '^) eV^ [IH]. To calculate the significances presented in the 
figure we have smeared the event distributions with oe = 0.2Eiy and ae = nip/ E„ as in 



Fig.ITSl 
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FIG. 20: Difference of significances in the (-E^-cos 0'') plane between the Am|^ uncertainty 



(Fig. 19) and the mass hierarchy (Fig. 13). 



V. CONCLUSIONS 

The main goal of our paper was to attract attention to the possibility of determination of 
the neutrino mass hierarchy with huge atmospheric neutrino detectors like PINGU, to pro- 
pose a method of quick estimation of the sensitivity of such detectors to the mass hierarchy 
and to outline challenges on the way of realization of this idea. 

1. After the determination of the leptonic 1-3 mixing angle, the structure of the neutrino 
oscillograms discussed in refs. pLl] and [TB] (see also Fig. [3]) is well determined, and the 
position of the main structures in the — cos 6z (kinematic) plane is fixed. For E,y > 1 
GeV these include the MSW resonance in the Earth's mantle domain as well as the MSW 
resonance and the three parametric ridges in the core domain of the oscillograms. 

2. The multi-megaton ice (water) detectors like PINGU will allow one to reconstruct 
the oscillograms and determine yet unknown neutrino parameters: the mass hierarchy (the 
sign of Amg^), the deviation of the 2-3 mixing from the maximal one, and in principle, the 
CP-violation phase as the next step. In addition, once the neutrino mass hierarchy has 
been established, the data from multi-megaton detectors should allow a significantly better 
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FIG. 21: Left panel- Significance of determination of the hierarchy as a function of the "fit" values 
of Am^g when the "true" value of Am^n^^^ = 2.35 • 10~^ eV^ (vertical line). Right panel - 
Significances as functions of Am|^ ^^^g after minimization as in the left panel (solid lines) and when 
Am23 = Am|;^ j^^g (dashed lines). A systematic uncertainty / = 5% and 5-yr PINGU data was 
used. 

determination of the value of Am|^. 

3. At the probability level the effects of the change of the mass hierarchy and of variations 
of 023 can be of order 1. However, there are several factors that substantially reduce the 
effects at the level of observable events. We have identified and studied in detail the following 
factors: 

a) u — u summation, which is related to the presence of both neutrinos and antineutrinos 
in the original neutrino flux. The hierarchy asymmetry survives due to a factor of ~2 
difference of the neutrino and antineutrino cross-sections as well as some difference of the 
original neutrino fluxes. 

b) Flavor screening, which leads to the suppression factors like (s|3r — 1) and (^23^ — 1) 
and is related to the presence of both Ue and z/^ in the original flux. 

c) Dilution of the significance: often large significances of the hierarchy asymmetry and 
other event number differences appear in bins where numbers of events are small. Then 
summation of signals in these high- significance bins and in bins where the numbers of events 
are larger but the significance is lower leads to a dilution of the high significance. 
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d) Smearing: finite energy and angular resolutions mean that integrations over rather 
large domains in the — 9^ plane have to be performed. These domains usually contain 
regions with different signs of the considered effect (e.g., hierarchy asymmetry). Therefore 
the integration leads to a significant suppression of the studied effects. 

e) Parameter degeneracies. 

4. We presented the significance plots for the determination of the neutrino mass hierar- 
chy, deviation of the 2-3 mixing from the maximal one and CP-phase. 

5. To evaluate the possibility to establish the neutrino mass hierarchy, we performed 
smearing of the event number distributions using Gaussian functions for reconstructing the 
true neutrino energies and zenith angles. Wc have studied the dependence of the significances 
integrated over certain ranges of and cos 9z on the widths of the reconstruction functions. 
Our preliminary estimates show that after 5 years of PINGU 20 operation the significance of 
the determination of the hierarchy can range from ~ 3a to lOcr (with parameter degeneracies 
taken into account) , depending on the accuracy of reconstruction of the neutrino energy and 
zenith angle. 

6. The smearing procedure implemented in this paper captures the main uncertainties of 
reconstruction of the true neutrino energies and zenith angles rather accurately. The smear- 
ing we have adopted gives a good estimate (at least at this stage of knowledge of future 
experimental characteristics) of the accuracy of reconstruction of the neutrino parameters. 
By varying the smearing parameters in rather wide ranges we covered essentially all possi- 
bilities of practical interest and explored the detecor resolutions necessary to achieve a given 
significance of the hierarchy determination. 

7. The parameter degeneracy effects can be significant, so that similar patterns of event 
distribution in the E^, — 0^ plane or in its parts can be obtained due to either changing 
the neutrino mass hierarchy or due to variations of the 2-3 mixing or of Am|]^ within their 
currently allowed ranges. For each parameter we identified the kinematic regions of the 
smallest degeneracy (where the effect we are interested in dominates). The effects of the 
parameter degeneracy can be reduced by selecting particular regions of integration over E^, 
and cos 9z. In addition, forthcoming measurements of neutrino parameters (in MINOS, T2K, 
NOvA and in reactor experiments) should provide us with more accurate values of these 
parameters and further reduce the effects of parameter degeneracy. 
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Appendix. Proof of the relation = Pa^(A — )• —A) in the limit Arn^i — )■ 

Let us prove that in the limit Am'^i — ?■ the oscillation probabilities for antineutrinos are 
given by those for neutrinos with the substitution Am|]^ — > —Aml^. A slightly different 
proof of this statement can be found in [29] . 

The evolution equation for the neutrino state vector in the propagation basis is 



i — i' = H{x)i> , 

ax 



(32) 



where v = {ug, z/g, ^3) and the effective Hamiltonian H{x) in the limit Amgi — )■ takes the 



form 



H(x) 



( s\^/\ + V{x) S13C13A ^ 




V 



S13C13A C?3A 



(33) 



^13^ / 

Here A = /\m\^/2E. The evolution matrix (the matrix of the transition amplitudes) then 
has the form 



A{x) 



( 



A,, A, 



e3 



1 



^3e ^33 



(34) 



From unitarity of this matrix it follows that \A-^^ = \A^r^\. For antineutrinos, one has to flip 



the sign of the potential V{x) in eq. (33). Obviously, this is equivalent to flipping the sign 
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of A in (33) and additionally changing the overall sign of the Hamiltonian H{x). The latter 
can be compensated by complex conjugating the evolution equation (32). Thus, we find 



A{u,A) = A*{u,-A). (35) 

Recall now that in the limit Amgi — ?■ the probabilities of various flavour transitions are 
expressed only through Pa = l^esP cos0x where 0x = ci^di^ii^l^) — (^^9{^l^) (see 
eqs. ([T])-([5])). Therefore we find that the oscillation probabilities for antineutrinos are given 
by those for neutrinos of the opposite neutrino mass hierarchy [i.e. with A — > — A). 
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